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Abstract 

We re-analyze the conditions for the phenomenon of intermittency (self- 
similar fluctuations) to occur in models of multifragmentation. Analyzing two 
different mechanisms, the bond-percolation and the ERW (Elattari, Richert 
and Wagner) statistical fragmentation models, we point out a common quasi- 
gaussian shape of the total multiplicity distribution in the critical range. The 
fixed-point property is also observed for the multiplicity of the second bin. 
Fluctuations are studied using scaled factorial cumulants instead of scaled 
factorial moments. The second-order cumulant displays the intermittency 
signal while higher order cumulants are equal to zero, revealing a large in- 
formation redundancy in scaled factorial moments. A practical criterion is 
proposed to identify the gaussian feature of light-fragment production, distin- 
guishing between a self-similarity mechanism (ERW) and the superposition 
of independent sources (percolation). 
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1. Two classes of multifragmentation mechanisms: percolation and ERW 
models 

In heavy-ion collisions, it is known that heavy excited systems are created and break into 
many lighter fragments of different sizes. This phenomenon is called multifragmentation. It 
occurs when the energy deposited in the system is sufficient to develop instabilities but not 
enough to totally evaporate the system into nucleons. A statistical analysis of the average 
spectrum is consistent with a power-law dependence of the fragment-size distribution. Some 
interesting interpretations of multifragmentation based on a phase transition in nuclear 
matter have been motivated by the experimental form of the fragment-size distribution 
I]!]!!, compatible with a power-law. Indeed, several critical phenomena, such as for instance 
percolation and liquid-gas phase transition exhibit a similar power-law of the fragment-size 
distribution at the critical point J3],§]. In these systems, the role of fragments is played 
by connected clusters of a given phase inside the other phase and the size distribution of 
clusters is proven to take the general form: 

N(s,£)~s- T -f(sE°) T,a>0, (1) 

where s is the cluster size and e characterizes the distance from the critical point. In a 
thermal phase transition e = T — T c with T c the critical temperature. In a bond-percolation 
model e = r c — r, with r the probability for the bond between two neighbouring sites to 
be broken and r c its critical valueQ. The genuine phase transition is well-known to only 
occur when the system is infinite. However finite systems show specific features related to 
the continuous limit, including finite-size corrections |||J. In multifragmenting systems, 
the parameter e is not well identified. It could be associated with the excitation energy 
deposited into the system by the collision. 



1 In multifragmentation problems, it is more convenient to use the variable r = 1 — q, where q is 
the usual bond probability j| . 
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One typical property of phase transitions is the existence of large, scale invariant fluctu- 
ations near the critical point || . In order to extract a signal from fluctuation patterns, one 
has proposed || in nuclear physics a method based on factorial moments. This method was 
originally designed |J for the rapidity spectra of ultra-relativistic multiparticle production 
but can be extended to other problems. Dividing the phase space into M bins of size 5, the 
factorial moment of order p can be definedF] by 

Fp(6) = EfL<n s >v ' (2) 

where n s is the number of "objects" in the s th bin and the average is taken over the set 

of events. In high-energy physics, "Objects" are particles distributed as a function of the 
rapidity variable. In nuclear physics, they are fragments distributed as a function of the 
size. 

In particle physics, the interest of factorial moments has been to deconvolute the statis- 
tical fluctuations due to the limited number of "objects" observed by event. It has revealed 
the existence of self-similar fluctuations, called intermittency, that is: 

F p (5)oc5- a r, (3) 

where a p is the so-called intermittency exponent. For event-by-event fluctuations of multi- 
fragmentation spectra ||, one considers M bins of size 5=4j along the fragment-size axis (A 
being the mass of the excited system). When M is varied, a behaviour like (3) was observed 
in the fragment-size distribution of ^Aung nucleus multifragmenting in a nuclear emulsion 
0. The signal was compared with that obtained from a site-bond percolation model in a 
network of size 6 3 , showing a similar behaviour at the critical point. The behaviour (3) was 
thus interpreted phase-transition signal. 

However the application of the method to the fragment-size distribution is less straight- 
forward than for rapidity spectra of particles. Indeed, the expression (2) is dominated by the 



2 We choose here the same definition as in Ref. fl5fl 
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first bin, i.e. the one containing the lightest fragments. In particular, no fragment of mass 
greater than A/p can contribute to F p . This introduces an important bias in the analysis, 
and has to be properly taken into account. 

On a less technical ground, several remarks have made the phase-transition interpreta- 
tions doubtful 0. In percolation models, the signal disappears when the size of the system 
goes to infinity. Moreover, F P (S) seems to strongly depend on the shape of the total mul- 
tiplicity distribution. In particular, when selecting events with the same multiplicity, the 
signal disappears. 

The interpretation is made more puzzling when considering a different class of multifrag- 
mentation mechanisms which leads to an intermittency behaviour without relying on a phase 
transition: the ERW model |J. Though there is no quantitative multifragmentation model 
using this mechanism, its simplicity has made it a fruitful toy model for understanding the 
mechanisms behind intermittency. In this case, the rule is to create fragments of different 
sizes with a Monte-Carlo procedure as follows. The average fragment-size distribution is 
constrained to be a power-law distribution P(s) ~ s~ r , r being a parameter of the model. 
We note A the total mass. An event is specified by a set of random numbers rji, 772, • • • , r\ n 
(with < r]i < 1), which defines a set of successive binary fragmentations (see Fig. 1). r]i 
defines a first fragmentation into two fragments of mass a\ and A — a\, where ai, an inert 
fragment present in the final event, is implicitely defined by 

Eai— 1 „— r y^ai „— r 

Ei^ < ' h " ES^' (4) 

The event is completed by repeating the operation with the mass A — a\ (or A — X)* =1 «i 
after k generations) until the total mass is exhausted. This model has many interesting 
properties 0, in particular, an intermittent signal can be seen for 1.8 < r < 2.0. 

Using the two different classes of mechanisms, percolation and ERW models we will 
re-examine the problem of interpreting intermittency by asking the following questions: 

(i) Is there a common scale-invariant mechanism behind the different realisations of 
intermittency? 



(ii) Are there statistical tools which could be more suitable than factorial moments for 
the study of multifragmentation spectra? 

(iii) Can we distinguish between different mechanisms? 

The plan of our study is as follows. In section 2, we show in both models the existence of 
fixed points, where the shape of the multiplicity is stable when the size of the system goes to 
infinity. This shape is quasi- gaussian for the total multiplicity and bin-size invariant but not 
gaussian when selecting the second bin. In section 3, we show that factorial cumulants ||11|| , 
instead of moments, are best suitable for the study of gaussian fluctuations of the multiplicity 
spectra. In the final section, we propose a criterion to distinguish between models, together 
with a physical interpretation of the different mechanisms in action. 

2. Quasi-gaussian fixed points 

Let us introduce well-known coefficients which characterize the shape of a statistical 
distribution. For a given variable m, these are defined by 

< (m- < m >) (i+2) > 

T< = — ~ (51) ( 5 ) 

< [m— < m >) 2 > y 2 ' 

In particular, we will consider the skewness (71) and sharpness (72) coefficients^. For a 
gaussian probability distribution 71 = and 72 = 3. The number 7! (resp. 72) describes 
the main assymetric (resp. symetric) variation with respect to the gaussian. Interestingly 
enough, we have found fixed points in the parameter space, where the values of 71, 72 
evaluated for the total multiplicity nearly reach gaussian values. We call them quasi-gaussian 
fixed points. 

We have performed a systematic analysis of the multiplicity distributions in the ERW 
(Fig. 2) and percolation (Fig. 3) models. Using ERW, we have first plotted the values 
taken by 71 and 72 for the total multiplicity distribution as a function of r (see Fig. 2-a) 
for different values of A (50 < A < 500). Each of the points reported on these figures is 
calculated with statistics of about 10 5 events. We note the existence of a value r ~ 1.8 where 



3 Here, 72 should not be confused with Campi's notation 
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all curves intersect. Interestingly enough, this value is in the range where the intermittency 
behaviour has been noticed ||. The values of 71 and 72 at the "fixed point" are close 
respectively to and 3 but nevertheless slightly different. The total multiplicity distribution 
is thus quasi-gaussian. We also note the existence of a minimum for 72 near r ~ 2 which 
confirms that the intermittent signal occurs when fluctuations are maximal^. This study 
has revealed the existence of a scale-invariant property of the total multiplicity distribution. 

In order to understand its connection with the intermittency signal, an analysis at fixed 
total size but varying the bin size appears useful. We thus studied in more detail the 
multiplicity distribution in the second bin ri2(5) (Fig. 2-b). The results are shown for 71 
and 72 as a function of r. The curves correspond to varying values of the bin size (4 < 5 < 64) 
with a fixed value of the system size A = 128. There is again a clear evidence for a "fixed 
point". We have also verified that this fixed point remains when the system size increases. 
We note that the corresponding value of r is similar but slightly different than that of Fig. 
2-a. Nevertheless, it remains in the region where the intermittent behaviour is noticed. The 
same analysis for the first bin distribution n\{8) is not relevant because the values of 71 and 
72 are strongly dominated by the distribution of fragments of mass one. In order to study 
the scale-invariance properties as a function of the bin size 5, it is thus necessary to avoid as 
much as possible the statistical bias due to the dominance of mass-one fragments. Moreover, 



there could be good experimental reasons to avoid the contribution of lightest fragments flO 
which can be contaminated by other processes than nuclear multifragmentation (e.g. pre- 
equilibrium, secondary emission, etc.). 

The same analysis was performed with the finite-size bond percolation model (Fig. 3). 
One considers a 3-dimensional cubic lattice with a probability r of breaking a bond. The 



4 When 72 > 3, the distribution is sharper than a gaussian, and the fluctuations are weaker. On 
the contrary, the fluctuations are larger when 72 < 3 (for instance 72 is equal to 1.8 for a uniform 
distribution). In our study, we found a minimum of 72 near 2.5. 
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critical value of the bond-breaking probability r c is around 0.76 (0.7512 in the continuous 
limit [|J). We have first plotted the values of 71 and 72 for the total multiplicity distribution 
as a function of r for various sizes of the system (4 3 < A < 10 3 ), see Fig. 3-a. Each point 
corresponds to the same statistics (10 5 events) as for the previous study. There is evidence 
for a fixed point in 71 near r c . The fixed-point behaviour is less clear for the sharpness 
parameter 72. The values of 71 and 72 remain respectively near and 3 confirming the 
existence of a quasi-gaussian shape in the critical region of the parameter. In the same 
region, intermittency has been observed ||. In the second bin, the analysis of 71 and 72 
(Fig. 3-b) clearly indicates the existence of a fixed point in the same parameter range. 

The strong similarities between the two models and the location of the fixed point in 
the same parameter region where intermittency occurs, calls for a deeper analysis of inter- 
mittency. Indeed, the quasi-gaussian features of the total multiplicity distribution require a 
specific treatment of intermittency which goes beyond the use of factorial moments. 

3. Factorial cumulants 

The existence of quasi-gaussian distributions leads us to introduce new statistical tools 
for multifragmentation, namely the factorial cumulants. In particle physics, they have been 



introduced |TT| in order to combine the elimination of the statistical noise using the factorial 
form with the well-known property of cumulants. Cumulants are a-priori able to disentan- 
gle genuine higher-order correlations from combinations of lower-order correlations which 
appear in factorial moments. In particular, gaussian fluctuations are governed by 2-body 
correlations implying the vanishing of cumulants of order greater than two. Let us apply 
this tool to nuclear multifragmentation taking into account the specific normalisation used 
in formula (2). We first introduce the generating function of unsealed factorial moments 

M 

G(A) = ]T<A->, (6) 

s=l 

where A is an arbitrary parameter and the average < ■ > is performed over the set of events. 
For instance, formula (2) can be identically written as 
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/ d p G(X) \ 

F p (5)~ ^" ' 



(7) 

A=l 



(^ A) ) P 

Using the logarithm of the generating function, the general expression of the scaled factorial 
cumulants can be written 

( dP\nG(X) \ 
\ d\P J 



KM 



For instance, the first cumulants read 

K 2 (§) = g£i_( <^K " 1) > " <n s > 2 ) 



(8) 



A=l 



K 3 (S) 



EZi<n s > 2 

gjli (<n s (n s - l)(n s - 2) > -3 <n s (n s - !)><n s > +2 <n s >' 

EZi<n s > 3 



v fx , E^i (<n.(n. - 1)K - 2)(n, - 3) > -4 <n,(n 4 - l)(n, - 2) ><n, >) 



E:u<n s >* 

Efli (-3 <n s (n s - 1)> 2 +12 <n s (n s - \)><n s > 2 -6 <n s > A ) 

Y.f =l <n s >* 



Z^ s =l \— <-> ^ n s yn s — i) ^ Tiiv/t s ^t s -i;^v/t S x' — u^iig^ j , . 



It is interesting to note that the choice of normalisation (2) implies that the scaled factorial 
cumulants cannot in general be expressed as combinations of scaled factorial moments. 
If however, the first bin dominates the evaluation of moments, one obtains the following 
approximate relations: 

K 2 (6)~F 2 (8)-1, 
K 3 (5)~F 3 (8)-3-F 2 (8) + 2, 

K 4 (S) ~ F 4 (S) - 4 • F 3 (S) - 3 • (F 2 (5)f + 12 • F 2 (6) - 6. (10) 

Note that the direct determination of cumulants using formula (8) could be useful to avoid 
the errors on the combinations of factorial moments. 

Let us apply our formalism to the ERW and percolation models. We have first computed 
the factorial moments using definition (1), see Fig. 4-a and 4-b. The result obtained with 
the ERW model is displayed in Fig 4-a with a total mass A = 128 using the critical value 
t = 1.8. The figure 4-b corresponds to the same observables for the bond-percolation model 
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with A = 6 3 and r ~ 0.75 = r c . The curves reproduce (for slightly different parameter 
values) the results of refs. || and |9[]. Note that, for completion, we have also reported the 
ERW factorial moments for the value r = 2.3 corresponding to the slope observed in the 
critical region of the bond-percolation model. The results exhibit a nearly linear increase as 
a function of the bin size which is typical of an intermittent behaviour. 

We have calculated the factorial cumulants corresponding to formula (9) for the ERW 
(Fig. 4-c) and percolation (Fig. 4-d) models. Cleary, only the second cumulant K 2 (S) 
is significantly different from 0. The higher order factorial cumulants are nearly zero for 
all bin size 5. Note that due to the smallness and the change of sign of cumulants, it is 
not convenient to plot their logarithm. Indeed, the cumulants obtained for both models 
correspond to strong cancellations between the large and positive scaled factorial moments. 
This confirms that there is a redundancy of information in factorial moments which is solved 
by using cumulants. 

The figures 4-c and 4-d show a clear indication of a rise of K 2 (5) especially for the per- 
colation model. When compared to Fig. 5, where scaled factorial cumulants are displayed 
for the bond-percolation model outside the critical region, the behaviour of K 2 (5) is re- 
vealing intermittent fluctuations. Indeed, an intermittent singularity in factorial moments 
should also appear in the second cumulant since it cannot be absorbed by the lower order 
contributions. For small values of K 2 (5) the first relation (10) gives 

K 2 (5) ~ log(F 2 (5)) ex -a p log(<5). (11) 

Such a behaviour is observed in Fig. 4. It is interesting to note that the dominance of 
second-order cumulants is valid for all values of the parameter r in the percolation model. 

Considering the abovementionned properties of scaled factorial cumulants for multifrag- 
mentation mechanisms, a series of comments are in order: 

(i) The scaled factorial moments of multifragmentation models seem to contain redundant 
information, since only the second cumulant is significantly different from zero. In the 
framework of the approximation (10), which is valid since the first bin is largely dominating, 
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it means that the scaled factorial moments can all be expressed as functions of 1(2(8). 

(ii) In the percolation model, the cumulant analysis indicates that the same property still 
holds outside the critical region, see Fig. 5. However, it is only around the critical region 
r ~ 0.75 that one observes a rise of 1(2(6) specific of intermittency. 

(iii) A qualitative difference exists between the behaviour of K 2 (6) for the ERW and the 
percolation models. They take very different values for the largest bin size (8 = A), i.e. the 
second cumulants of the total multiplicity distributions are significantly different. 

(iv) We also remarked that the slope of 1(2(6) decreases with increasing system size, 
similarly to what has been observed for scaled factorial moments ||. However, the order- 
of-magnitude difference still remains between the values of K 2 (8) for the ERW and the 
percolation models. 

4. Discussion and summary of results 

Our analysis of multifragmentation spectra using the comparison of two generic mecha- 
nisms, ERW and percolation, calls for a physical interpretation. The existence of a scaling 
behaviour for the profile parameters 71, 72, being a common feature of both models in their 
respective critical regime, points to the existence of a scale-invariant property of multiplicity 
distributions. It is phenomenologically clear that this "fixed-point" is related to the inter- 
mittency signal and thus to a certain type of criticality of the system. The criticality is 
explicit for percolation, since it corresponds to a phase transition in the continuous limit, 
but is hidden in the formulation of the ERW mechanism. 

Observing at the same time a quasi-gaussian shape of the total multiplicity and a large 
dominance of K 2 (5) over higher order cumulants unravels a pronounced gaussian feature of 
multifragmentation mechanisms when small-mass fragments dominate the observables. Two 
main interpretations may naturally explain these features. It may come either from the 
existence of a self-similar fragmentation mechanism with a scale invariant gaussian fixed- 
point or from the superposition of many independent sources of fragments leading to a 
gaussian distribution through the central-limit theorem. As we once more want to stress, 
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the gaussian property is mainly inherent to the production of lighter elements. 

In order to distinguish between the two options, let us study the two parameters which 
define a gaussian distribution namely the distribution average < m > and its variance 
a =< (m— < m >) 2 >. One expects [[12] a j < n >^ cste for self-similar models and 
a I <n> oc l/y/<n> for the independent sources. Indeed, in this latter case the number 
of independent sources is expected to grow with the multiplicity. Then the central-limit 
theorem gives the quoted prediction for aj < n >. In Fig. 6, we display the quantity 
oj < n > for each model for both the total multiplicity distribution and the mass-one 
fragment distribution. For the total multiplicity, power-law fits give a j <n> oc <n > -001 for 
the ERW model and o j <n> oc <n>~ 0AA for the percolation model. Interestingly enough, 
the fits are very similar for the mass-one fragments with a power dependance <n>~ 0A9 for 
the percolation model very suggestive of a central-limit signature. Hence, the ERW is a self- 
similar mechanism, while percolation is gaussian with a central-limit type for the production 
of light fragments. 

Let us summarize our main results by proposing the following new tools of analysis for 
multifragmentation fragment distributions. 

(i) The shape parameters 71, 72 for the total multiplicity and for the second-bin multi- 
plicity can be studied as a function of the system size, and for the second function 
of the bin size. We expect a fixed-point behaviour for the critical region of multifragmenta- 
tion. Intermittency is expected to occur in the same region. Note that the second bin could 
be experimentally easier to study, since it avoids the ambiguity on the origin of lightest 
fragments and can be used at fixed total-size by varying only the bin-size. 

(ii) The study of fluctuations using scaled factorial moments suffers from a large redun- 
dancy of information. For the class of models under study, we suggest using instead the 
scaled factorial cumulants of fluctuations. 

(iii) The gaussian feature of the light-fragment distributions revealed both by the quasi- 
gaussian fixed point property and the scaled factorial cumulants is an important character- 
istics to be studied in nuclear multifragmentation. Note that this gaussian feature has also 

11 



been noticed in a different model fll3| , confirming the interest of studying this property in 
detail. Two mechanisms seem to give an alternative for explaining the gaussian features, 
namely either a self-similar fixed-point or a superposition of independent emission sources 
for light fragments. They can be distinguished by dependence of a j <n> with < n > or 
equivalently by the factorial cumulants for large bins. As an application, we found that the 
ERW model is self-similar while percolation probably leads to independent sources for the 
production of light fragments even outside the critical region. We expect the self-similar 



property to be shared by models based on a tree structure of multifragmentation ||14| . |15 



In the same spirit we also expect models based on second-order phase transitions to fol- 
low the same trend governed by the central-limit theorem as the percolation model for the 
distribution of light fragments. 
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FIGURES 




FIG. 1. Event generation in the ERW model. The set of successive fragments 01,02 ■■■ ,a n is 
denned by the set of random number 771, 772, • • • , r? n (see formula (4)). The generation of random 
numbers is stopped when Yh=i a % > A. 
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FIG. 2. Shape parameters 71, 72 of the multiplicity distributions for the ERW model. 
a) Top: 71, 72 (total multiplicity distribution) as a function of r for different system sizes 
(A = 50, 100, 200, 300, 400 and 500). b) Bottom: 71, 72 (second bin) as a function of r for different 
values of the bin size 5 (A = 128). 
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FIG. 3. Shape parameters 71, 72 for multiplicity distributions for the bond-percolation model. 
a) Top: 71, 72 (total multiplicity distribution) as a function of r for different system sizes A. b) 
Bottom: 71, 72 (second bin) as a function of r for different values of the bin size 5 (A = 6 3 ). 
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FIG. 4. Scaled factorial moments F p {5) and cumulants K p {8). a) \og(F p {$)) as a function 
of -log(<5) in the ERW model for A = 128, r = 1.8 (solid line) and r = 2.3 (dashed line), b) 
log(F p (<5)) as a function of — log(<5) in the percolation model for A = 6 3 and r = 0.75. c) K p {8) 
vs - log(<5) in the ERW model for A = 128 and r = 1.8 (solid line) and r = 2.3 (dashed line), d) 
K p (5) vs — log((5) in the percolation model for A = 6 3 and r = 0.75. 



18 



a) 



b) 



0.015 



0.010 



0.005 - 



0.000 



-0.005 




-log(S) 



0.002 p 

0.001 
0.000 
-0.001 - 
f-^- -0.002 - 
^ -0.003 - 



-0.004 



-0.005 - 



-0.006 



-5.5 



P=2 
p=3 
p=4 



• • • » — » > 



-3.5 



-1.5 



-log(S) 



FIG. 5. Scaled factorial cumulants K p (5) outside the critical region for the bond-percolation 
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FIG. 6. The scaled variance a/ < n> as a function of < n>. a j < n > is displayed as a 
function of the average multiplicity. Circles correpond to the variance and average of the total 
multiplicity distribution; stars correspond to the variance and average multiplicity distribution of 
mass-one fragments. For the ERW model with r = 1.8, each point corresponds to a given total 
mass A taken in the intervalle [100,10000]. Power-law fits are displayed by continuous lines, the 
powers are —0.01 (circles) and —0.06 (stars) in the ERW case. For the bond-percolation model, 
r = 0.76 and each point corresponds to a fixed value of A (A = 3 3 , 4 3 , • • • , 20 3 ). The power-law fits 
are found respectively equal to —0.44 (circles) and —0.49 (stars). 
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